Crystalline order in superfluid 3 He films 
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We predict an inhomogeneous phase of superfluid 3 He films in which translational symmetry is 
spontaneously broken in the plane of the film. This phase is energetically favored over a range of 
film thicknesses, D C2 (T) < D < D C1 (T), separating distinct homogeneous superfluid phases. The 
instability at the critical film thickness, D C2 « 9£(T), is a single-mode instability generating striped 
phase order in the film. Numerical calculations of the order parameter and free energy indicate 
a second-order instability to a periodic lattice of degenerate B-like phases separated by domain 
walls at D C1 w 12 £(T). The striped phase should be identifiable in transport and nuclear magnetic 
resonance experiments. 

PACS numbers: 67.57.-z,67.57.Bc,67.57.Np 



The phases of superfluid 3 He provide a beautiful 
example of spontaneously broken symmetry in con- 
densed matter physics, exhibiting properties common 
to superconductors, nematic liquid crystals and anti- 
ferromagnets. Many of the unique physical properties 
of superfluid 3 He, including the spectrum of low-energy 
excitations, are connected to the spontaneous breaking 
of orbital and spin rotation symmetries in combination 
with global gauge symmetry that is associated with su- 
perfluidity and superconductivity. Inspite of the com- 
plex order that develops, the bulk A and B phases of 
3 He are translationally invariant. Indeed translational 
symmetry is generally assumed to hold even in reduced 
dimensions, e.g. superfluid filmsi^ 

NMR measurements on relatively thick (/im) films show 
evidence of an A- to B-like transition predicted within 
the context of Ginzburg-Landau (GL) theory. 3 How- 
ever, unexplained anomalies in film flowi and third 
sound experiments 5 suggest that our current theoret- 
ical understanding of the phases of superfluid 3 Hc films 
is insufficient. One of the intriguing questions raised by 
these experiments is whether or not there may be qual- 
itatively new phases stabilized in reduced dimensions.- 

Here we report the theoretical prediction of a phase of 
superfluid 3 He exhibiting spontaneously broken trans- 
lational symmetry, i.e. crystalline order. This phase is 
shown theoretically to be the stable ground state of a 
superfluid 3 He film, with the broken translational sym- 
metry occurring in the plane of the film. The mecha- 
nism responsible for this phase is competition between 
surface depairing and domain wall formation between 
degenerate ground states, and is generic to 3 He con- 
fined in at least one spatial dimension. 

The superfluid phases of 3 He are Bardeen-Cooper- 
Schricffcr (BCS) condensates of orbital p-wave [L = 1) 
Cooper pairs formed from quasiparticles with zero total 
momentum (+p, — p) near the Fermi surface in spin- 
triplet (5=1) states. In terms of the basis of triplet 
states the order parameter is given by 

A = A+(p)| Tt)+A-(p)| U) + Ao(p)-^=| TI + IT), 

(1) 



where A m (p) = J2i= x , y , z A mi Pi f° r m = 0,±1. There 
are two bulk phases of superfluid 3 He in zero field. 
For a narrow temperature range near T c at high pres- 
sures, p > p c — 21 bar, 3 He condenses into the A phase 
with an order parameter of the form, A + = A_ = 
and Ao = A(T) (p x + ip y ). This phase exhibits anti- 
ferromagnetic spin correlations, and an orbital state 
that breaks time-inversion symmetry, i.e. a conden- 
sate of pairs with orbital angular momentum +h. The 
B-phase, which is the stable state over most of the 
phase diagram in zero magnetic field, is a superposi- 
tion of all three triplet spin states and all three or- 
bital states, with A + = A(T) (p x - ip y ) /^2, A_ = 
A(T) (pa; + ip y ) /V2, A = A(T)p z . This state de- 
scribes a condensate of spin-triplet, p-wave pairs in a 
state with total angular momentum J = 0. There is a 
continuous manifold of B-phase states related by a rel- 
ative rotation of the spin and orbital coordinate axes. 
Surface and nuclear dipolar energies resolve most, but 
not all, of the degeneracy. In addition to the bulk A- and 
B-phases, the planar (P) phase is a possible ground state 
for thin films of 3 He. The P-phase is a two-dimensional 
version of the B-phase with Ao = 0. Alternatively, the 
P-phase is an equal amplitude superposition of degener- 
ate, time-reversed A-phase orbital states with opposite 
angular momenta. As a result the P-phase is degener- 
ate with the A-phase in the weak-coupling BCS theory, 
but preserves time-inversion symmetry. 

Here we consider 3 He films of uniform thickness, D, 
bound to a solid substrate. The liquid-vapor inter- 
face is assumed to be perfectly reflecting and atomi- 
cally smooth. Thus, we consider p — > bar. This is 
also the weak-coupling limit for superfluid 3 He, as in- 
dicated by the heat capacity jump AC/Cn — * 1-43 for 
p — ► Obar. 7 Substrates may provide a range of scat- 
tering from specular to diffuse scattering depending on 
the degree of roughness. We consider specular and fully 
diffuse scattering using boundary conditions described 
in Ref. H 

The order parameter for the film geometry is defined in 
terms of x- and y-axes which lie in the plane of the film 
and the z-axis perpendicular to the film. Scattering of 
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FIG. 1: Left panel: Q z vs. T phase diagram showing reen- 

trance for homogeneous phases. Right panel: Instability on- 
sets for higher Q z (T) [orange line with boxes] for an inhomo- 
geneous phase with in-plane wavevector Q X (T) [red dashed 
curve]. Note that £o = Hvf /^ttTc — 77 nm at p — bar. 

3 Hc quasiparticles off the free surface and substrate sup- 
presses the orbital Pz-component of the order parame- 
ter (for either specular or diffuse scattering) in films less 
than about 1/im thick. In such thin films the A-phase 
or the P-phase is stable; in the weak-coupling limit 
these states are degenerate even with strong pairbreak- 
ing from diffuse scattering. In thin films the excitation 
spectrum is typically dominated by gapless excitations. 
The effect of this spectrum on strong-coupling energies, 
combined with relative importance of gradient energies 
for low-temperature, low-pressure thin films means that 
we cannot infer the relative stability of phases from 
what is known about strong-coupling energies in bulk 
3 He. 

Microscopic calculations show that as the film thick- 
ness increases surface scattering is unable to com- 
pletely suppress the p z -component of the order param- 



eter. Results by several groups predict that equilib- 
rium phase for film thickness, D > D C (T) ~ 10£(T), 
is the deformed B phase described by the order pa- 
rameter, A B = (A||p x , A||p y , A z p z ) f iS with A z = 
Aj_(T) sin(7r z/ D)m&&The transition is first- or second- 
order at a critical film thickness, D C (T), depending on 
whether the low temperature phase is the A-phase or 
the P-phase. The phase boundary, taken from our ear- 
lier calculation^ is shown in the left panel of Fig. ^ 
in terms of the critical wavevector, Q Z (T) = w/D c (T). 
The inset emphasizes the re-entrance (A <— B <— A) for 
T < 0.42 T c near the critical line, which suggests that a 
lower energy state at low temperatures, in the vicinity 
of the critical line, may be achieved by an inhomoge- 
neous phase that incorporates features of both phases. 
This is the case, but as we show below the structure of 
inhomogeneous phase is more complex than any of the 
homogeneous phases and evolves over a relatively wide 
range of film thickness. 

The transition from B-phase to the P-phase is second- 
order on the critical line, D C (T). Thus, we first look for 
a second-order instability to an inhomogeneous phase 
that pre-empts the P-B transition. Our starting point 
is the weak-coupling gap equation, 

HI) a: - r " l^-- 

xT£ ffW(p-,R;e m )-, A(± . )( ^ R) l , (2) 

m V |£m| / 

where A' •'(p, R) are the real (+) and imaginary (— ) 
parts of the order parameter, and f( '(p, R;e m ) = 
[f(p, R;e TO ) ± f (p, R; — e m )*]/2 are the corresponding 
pair propagators in the Matusbara formulation for 
equilibrium Fermi superfluids^i^ These objects satisfy 
second-order mode equations with the order parameter 
providing source terms, 



I( V/ .V) 2 f (+) -^f (+) = -™ m fA< + » + j' + )l+- [ A ,+) (A (+ »^ l+) ) + A<-»(A( + ».^) 
4 L J oj m 

- A M x (A ( ~> x 8 {+) ) + A (+) x (A ( ~' x < 5 ( " ) )] , (3) 
i(v / -V) 2 f ( - ) -^f"=- 7 ro, m [a<->+«<->1+-^ [ A(-'(A(-».«H ) + A (+) (A (-).#)) 

4 L J LU m 

- A (+) x (A (+) x 8^) + A ( - } x (A (+) x <5 (+) )] . (4) 



Note that v/ = u/p is the Fermi velocity, e m — 
(2m + l)7rT is the Matsubara energy, and uj^ = e 2 m + 
|^(+)|2 _|_ |^(-)|2 j wnere ^y(±) i s the order parame- 
ter of the unperturbed, translationally invariant phase. 
Lastly, o^' is the first-order correction we seek to find. 
These equations are valid up to first order in order pa- 
rameter corrections, and in their derivation we assumed 



that the gradient terms of in strongly confined 
space are of the same order as f themselves. 

For the P-state we can fix the overall phase so that 
A is real. We then have = and A (+) = 

A||(z) (PzjPyjO). The instability to an inhomogeneous 
phase is then a single-mode instability for pairs with 
zero spin projection along The eigenfunction for 
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the instability has the form, 

*,(p,R)= J2 a ^(Q) eiQ ' R Pj ' (5) 

j=x,y,z 

For a single mode instability in the plane of the film we 
can choose Q = (Q x , 0, Q z ). The resulting solution for 
the Fourier component of the S z = pair propagator 
is, 

/i+)(p,Q) = ^(Wy^^- • (6) 
l( V/ -Q) 2 +<(p) 

Application of the boundary conditions to the pair 
propagator at the free surface (z = D) and substrate 
(z = 0) yield a set of eigenvector equations for the un- 
stable mode. Here we discuss only specular scattering 
by the substrate, in which case the boundary condition 



appearance of the order parameter ampli- 
tudes, a z ,x{Qx, Qz) exp(iQ x x) cos(Q z z) and 
o-z,z{Qx, Qz) e~xjp(iQ x x) sm(Q z z). These amplitudes are 
shown in the left panel of Fig. [21 for T — 0.5 T c and 
D = 9.3 £o < D C (T). Note that the full solution for 
the order parameter above the lower critical thickness 
also shows very small oscillatory amplitudes for the 
in-plane spin-components, e.g. a XjZ . 




0.1 








0.05 





a xz 


/ a zz 


0.05 






— - "if 


-0.1 









-10 



f (± \x,z = 0,p) = f (±) (a 



0,P) 



(7) 



for any x and p, with p = p — 2n(n • 
p); and similarly for the free surface. The 
boundary conditions reduce to connections between 
the Fourier components of the order parameter, 

Q>z,x(,Qxi Q z) ~\~Q>z t x\Qxi Qz) ; Q>z,y\Qxi Qz) 

+a z ,y(Q x ,Qz) ,a z ,z(Qx,-Qz) = -a ZtZ (Q x ,Qz) , and 
fixes the wavevector Q z = n/D in terms of the film 
thickness D at the instability. These results and the 
gap equation generate the eigenvalue equations for the 
mode amplitudes, a z ,j(Q), 



ln(T/T c ) a z 



E 

j=x,y,z 



, i = x,y,z, (8) 



with the matrix elements given by 



47T 



l(v/ 



(9) 

Translational symmetry is unbroken along the y-axis 
in which case I xy = I yz = 0. The mode ampli- 
tudes separate into linearly independent blocks: a 2D 
(a z>x , a>z,z) an d a ID (a z . y ) block. A non-trivial so- 
lution to Eq. exists if, (ln(T/T c ) - I yy ) = or 
{{\n{T/T c )-I xx ){\n{T/T c )-I zz )~I 2 xz } = 0. The eigen- 
value equation for the ID mode amplitude has a maxi- 
mum unstable wavevector only for the transition to the 
homogeneous phase, Q z = tt/D c (T), Q x = 0. How- 
ever, the eigenvalue equation for the 2D block gives an 
unstable mode Q Z {Q X ,T) that pre-empts the homoge- 
neous transition. The maximum value of Q Z {Q X ,T) as 
a function of Q x for each temperature determines the 
lower critical film thickness, D C 2(T) < D C (T), for the 
transition to an inhomogeneous film with broken trans- 
lational symmetry in the plane of the film. The critical 
wavevector, Q Z (T), and the locus of values of Q X (T) 
are shown in the right panel of Fig.Q. 

The key signature of spontaneously broken 
translational symmetry the xy-plane is the 



FIG. 2: Order parameter amplitudes (in units of 2ttT c ) for 
the stripe phase at T = 0.5T C along the film for z w 2.5£o. 
Left panel: D c2 < D = 9.3 < D c . Right panel: D c < 
D = 10£o<-Dd. 

Modes with the same magnitude for the unstable 
wavevector, Q x , but different orientation in the plane 
of the film are degenerate. In the absence of an exter- 
nal bias to select the direction of the unstable mode, 
the instability may propagate in any direction in the 
plane of the film. For D > D c2 the spatial structure 
of the order parameter that is realized is determined 
from the minimum free energy. This phase may ex- 
hibit one-dimensional, stripe-phase order, or possibly a 
two-dimensional structure defined by two non-collinear 
wavevectors, e.g. a triangular lattice. A comparison 
of the possible minimum energy configurations of the 
inhomogeneous phase has not been carried out. Here 
we focus on the structure of the one-dimensional stripe 
phase. The broken symmetry phase persists for film 
thickness well above the original critical line, D C (T), 
for the homogeneous A-B transition. For thicker films, 
actually for D > _D c2 (T), the gap equation includes 
nonlinear driving terms that couple modes with differ- 
ent wavevectors. The ground state is periodic, but the 
structure is non-sinusoidal. The right panel of Fig. [21 
shows the order parameter amplitudes for a film with 
D = 10 £o > D C {T). The basic structure of this phase is 
indicated by the amplitude a z . Zl which has developed a 
soliton-like structure separating "domains" of degener- 
ate B-like phases: e.g. = (Ayp^, Aypy, — Aip z ) and 
Aj = (Anpas, A||p w ,+Axp s ). Also, centered on the soli- 
ton is a non-B-like phase, represented by a ZjX , bound to 
the domain wall. 

This basic structure also provides a clue to the under- 
lying mechanism stabilizing the inhomogeneous phase; 
it is the competition between the energy associated 
with surface pairbreaking and the energy cost of a do- 
main wall separating two degenerate B-like phases>ii 
Consider the two trajectories (labelled 1 and 2) shown 
in Fig. [21 The left panel shows a homogenous B- 
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phase, while the right panel shows two degenerate B- 
like phases corresponding to amplitudes — Aj_ left of a 
domain- wall and +Aj_ to the right. 




FIG. 3: Left panel: Surface reflection (trajectory f) with 
p z — > — p z leads to strong pairbreaking for the homogenous 
B-phase. Right panel: Trajectory 2 is a the strong pair- 
breaking trajectory because of the sign change Ax — > -Ai. 

For the trajectory 1 that reflects from the free surface 
we have p z — > — p z . This sign change is the origin of 
pairbreaking by a specular surface in 3 He-B; it leads the 
suppression of Aj_ and the formation of surface Andreev 
bound states^ The energy cost is directly related to the 
spectrum of surface states. By contrast the shallow tra- 
jectory which passes through the dashed plane without 
intersecting a surface encounters a nearly uniform order 
parameter resulting in little or no pairbreaking. 

Surface pairbreaking can be suppressed locally by com- 
pensating the sign change that results from surface re- 
flection. In particular, for the domain wall configura- 
tion shown in the right panel of Fig. [21 trie sign change 
for the scattering trajectory (p z — > — p z ) is compensated 
by the sign change associated with the degenerate states 
on opposite sides of the domain wall (±A Z ). Thus, over 
a few coherence lengths near the domain wall, surface 
reflection does not lead to strong pairbreaking, and cor- 
respondingly the energy cost of surface scattering is re- 
duced. However, it is not all "savings". There is an 
energy cost for the domain wall. The shallow trajec- 
tory crossing the domain wall now incurs a sign change. 
Pairbreaking occurs near the domain wall and a spec- 
trum of Andreev bound states forms on the interface. 

For very thick films (D ^> D C (T)) the translationally 
invariant B-phase is favored because the surface pair- 
breaking energy is small compared with the pairbreak- 
ing cost of a domain wall. But, for sufficiently thin 



films the domain wall energy is less than the surface 
pairbreaking energy and the broken symmetry phase is 
favored. The critical line where one domain wall is fa- 
vored over the uniform B-phase is D c i(T). For D < D c \ 
multiple domains are favored. Further reduction in the 
film thickness favors more domain walls until they dis- 
solve into the P-phase at the D C 2(T), or a first-order 
transition to the A-phase occurs. 

Figure 0| shows the T vs. D phase diagram for 3 He films 
at p — bar. Calculations of the critical film thickness 
D c i (T) based on the formalism described in Ref. |8| give 
a range of film thicknesses of order 0.75 /im < D < 
1.0 pm for T < 0.75 T c . 

In conclusion, our calculations predict that films of 
superfluid 3 He should exhibit an inhomogenous phase 
with spontaneously broken translational symmetry in 

1.0 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 




FIG. 4: T — D-phase diagram for superfluid 3 He films in the 
weak-coupling limit (p = Obar). The phase with crystalline 
order separates translationally invariant A(P) and B phases. 



the plane of the film over a substantial range of temper- 
atures and film thicknesses. This phase has no analog in 
bulk 3 He, and should be identifiable by its anisotropic 
transport properties. For example, the in-plane thermal 
conductivity should exhibit a reduced heat conductiv- 
ity normal to the direction of the stripes. Signatures of 
the inhomogenous phase should also be observable as a 
broadening of the NMR linewidth. 
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